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Abstract. We present recent BABAR results on searches for dark photons, long-lived scalar 
particles and new 7r°-like particles. 


1. Introduction 

In 2010, the Pamela experiment reported a positron excess for energies greater 10 GeV that 
is increasing with energy [I]. The excess was confirmed by the Fermi Large Area Telescope in 
2012 |2|. Two years later, the AMS experiment published higher-precision data up to 300 GeV 
confirming an excess of the positron fraction in the 10 - 250 GeV energy range |3|. Figure [I] 
shows the measurements from Pamela, Fermi and AMS. In the same energy range, Pamela and 
AMS did not observe any antiproton excess. Though an astrophysical explanation is very likely, 
theorists have come up with new dark-matter scenarios favoring light particles. In the past, 
dark-nratter particles have been associated with new heavy particles predicted in extensions of 
the Standard Model (SM), such as neutralinos in supersymmetric models. Physicists have been 
searching for weakly-interacting massive particles (WIMPs) in dedicated experiments for many 
years [4 - 201. The observation of the positron excess has triggered searches for light dark matter 


in different accelerator experiments. 



Figure 1. Positron excess observed by Pamela, Fermi and AMS. 

The SM may be connected to the dark sector through so-called portals [2l]. These are lowest- 
dimensional operators that may provide coupling of the dark sector to the Shfl Table 0 lists 
these operators for couplings to new vector particles, pseudoscalars, scalars, and neutrinos. At 
low-energy e + e _ collisions, the light vector portal is the most accessible portal, but the scalar 


1 higher-dimensional operators are mass suppressed. 
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portal can be probed as well 
symmetry is introduced 


25 . In the simplest realizations of the dark sector, a new U{1)' 


26 . The associated gauge boson, called dark photon (. A '), couples to the 


weak hyper charge with a mixing strength e. An effective interaction between the dark photon 
and the electromagnetic current arises after electroweak symmetry breaking. For an Abelian 
interaction, this mechanism would explain WIMP annihilation into SM fermions 125 27,28 


To 

accommodate the recent anomalies in cosmic rays, the dark photon mass is constrained to lie in 
the MeV/c 2 to GeV/c 2 range. 

In recent years, we have searched for new physics in e + e _ collisions collected with the BABAR 
experiment at center-of-mass (CM) energies around 10 GeV 
Higgs h! plus a dark photon A' in the process e + e~ —> h!A! 


29, 30 . We searched for a dark 
A'A'. The cross section is 


h' 

proportional to the product of fine structure constant ap and e 2 . Assuming a flat prior, we set 
90% confidence level (CL) Bayesian upper limits on an -e 2 of 10 ~ 10 to 10 -8 for h! and A 1 masses 

We also searched for 
ji ~, fi + fi~y + fi~ [321. 


in the 0.8 - 10 GeV/c 2 and 0.25 - 3 GeV/c 2 regions, respectively [31 
dark non-Abelian gauge bosons in e + e~ —> W'W' —> e + e~e + e ~, e + e~ y 

Furthermore, we searched for a new pseudoscalar A 0 in Y(2 S,3S) radiative decays: A { 

For the y + y 
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t + W [35], hadrons [36] , and invisible particles |37j. For the y + y~ decay mode, 
we set 90% CL Bayesian upper limits on the branching fraction of 0.26 x 10 ~ 6 to 8.3 x 10 ~ 6 
for A 0 masses between 0.21 GeV/c 2 and 10.1 GeV/c 2 33 . For the other modes, we obtained 

gg , ss setting 90% CL upper 
for A 0 masses in the 1.5 - 9.0 GeV/c 2 


similar results. We also studied the decay Y(1S) —> 7 A u , A 0 
limits on the branching fraction at a level of 10~ 6 to 10” 2 
range [38]. 


Table 1. Portals for connecting the dark sector to the Standard Model where F^ v is the 
strength of the Yang-Mills field, F lflu is the strength of the dark photon held, e is the mixing 
parameter of a dark photon to a real photon, a is the axion held, f a is the scale at which 
the Peccei-Quinn global 17(1) symmetry is spontaneously broken 122 


24 , G is the gluon held 


strength, S is a dark scalar held with coupling strengths (A and y) to the Higgs held H , and N 
is the sterile neutrino held with coupling to the Higgs held. 


Portal 

Particle 

Operator 

” Vector” 

” Axion” 

” Higgs” 

’’Neutrino” 

Dark photons 

Pseudoscalars 

Dark scalars 

Sterile neutrino 

eF^F^ 

{llS + A S 2 )H^H 

VnLHN 


2. Search for Dark Photons 

In the mass region below 10 GeV/c 2 , e + e - collisions provide an excellent laboratory for 
producing dark photons A' [39j[40|. The processes e + e _ —> 'yA' —> r yl + (.~ (jA = or gA) 
depicted in Fig. [ 2 ] are well suited since they have high sensitivity for A 1 production because of 
high selection efficiencies and large A' decay rates as Fig. [ 3 ] shows [39] . The cross section with 
respect to e + e _ —> 77 is reduced by e 2 , where the mixing parameter e is expected to lie in the 
range of 10 -5 to 10 ~ 2 . 

We have studied e + e _ 7 and n + 7 hnal states at CM energies larger than 200 MeV 
using an integrated luminosity of 514 fb _1 of e + e _ collisions 41 . The data were recorded 
at the T(45), T(35), T(25) peaks and 40 MeV/c 2 below the T(4S') peak with the BABAR 
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Figure 2. Lowest order diagram 
for e + e _ —»• 7 A' —> ^l + l~ ^qq. 



Figure 3. Branching fraction predictions for A' 
i + l~ and A' -7 qq as a function of uia'- 


detector [29, 30 at the PEP-II asymmetric storage ring at the SLAC National Laboratory. 


We select two oppositely-charged leptons plus a photon and apply a constrained fit to the beam 
energy and interaction point (IP). We use additional kinematic constraints to improve purity, 
require good quality on the photon, electron and muon and remove electron conversions. We 
simulate the background processes e + e“ —> e + e~('y ) and e + e _ —> 77 ( 7 ) with the generator 
BHWIDE 42 , e + e _ —> with the generator KK2F [431 and resonance production in 

initial state radiation e + e _ —> 'yX (X = J/0, '0(2S’), Y(IS'), Y(25)) using structure function 
techniques [44[|45|. For the simulation of the detector response, we use GEANT4 


46 


Figure |4| (left) shows the dielectron invariant-mass spectrum in data and in simulation 
with the generator BHWIDE. Generally, the agreement between data and Monte Carlo (MC) 
simulation is good except for the low-mass region where the BHWIDE generator fails to 
reproduce the m ee invariant-mass spectrum. On the other hand, the MADGRAPH generator [47 
reproduces the low m ee region in a consistent way but the limited sample size introduces too 
large uncertainties. Therefore, our signal extraction does not rely on the MC simulation. 



Figure 4. Top: measured m ee (left) and tur (right) invariant-mass spectra in comparison to 
simulations; bottom left and right: the data/MC ratio. 


For the dimuon sample, we define the reduced mass tiir = — 4m^) since the turn-on 

near threshold is smoother. Figure [4] (right) shows the vir distribution in data and simulations. 
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The KK2F generator |43| reproduces m,R spectrum reasonably well. For both modes, we 
determine the selection efficiency from simulation yielding e ee7 = 15% and e MAt7 = 35%. The 
efficiency reduction in the eey mode is due to mainly pre-scaling of radiative Bhabha events in 
the trigger. However, radiative Bhabhas are still the dominant background. Thus, the sensitivity 
for the dark photon search is dominated by the fj, + /x - mode. We extract the signal yield as a 
function of niA' by performing independent fits to the m ee and itir mass distributions for each 
beam energy covering the range 0.02 < tiia 1 < 10.2 GeV/c 2 and 0.212 < uia' < 10.2 GeV/c 2 
for the eey and ///ry mode, respectively. In these fits, we vary the step size to about half the 
value of the expected mass resolution a mA , and cover a mass range that is at least 20a mA ,. The 
signal resolution is estimated by fits to a Gaussian function for various A! mass values and is 
interpolated to all other masses. The typical steps sizes vary between 1.5 MeV/c 2 and 8 MeV/c 2 
yielding 5704 (5370) mass hypotheses for the eey (yu/ry) channel. 



Figure 5. Measurement of a(e + e — > yW — y 7 e + e ) (left) and significance distribution (right). 
Grey vertical lines show the excluded regions around the vector mesons. 



m A (GeV) -4 -3-2-10 1 2 3 4 

b s 

Figure 6. Measurement of a{e + e~ —> yvT —> 7 (left) and significance distribution (right). 
Grey vertical lines show the excluded regions around the vector mesons. 

The likelihood function used to fit the observed spectra contains a signal com¬ 
ponent, radiative dilepton background and peaking background from vector mesons 
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(w, (j), J/V’, ip(2S), Y(1 S), T(2S)). The signal probability density function (PDF) is modeled 
directly from the simulated signal mass distribution by a non-parametric kernel pdf and is inter¬ 
polated between the known simulated masses using an algorithm based on the cumulative mass 
distributions 48 . We estimate a systematic uncertainty of 5% - 10% for this procedure. The 
PDFs for the radiative dilepton backgrounds consist of higher-order polynomials. Though we 
parameterize the vector mesons ui, (j), J/4>, ip(2S), T(IS') and T(2 S) with Gaussian and Crystal 
Ball functions, we exclude the mass regions around them in the dark photon search. We de¬ 
termine the e + e _ —> 7 A' —> r ye + e~ and e + e _ —> 7 A 1 —> 7 n + cross sections by dividing the 

extracted signal yield by efficiency, luminosity and dark photon decay branching fractions 1 39 1 . 
Figures [ 5 ] (left) an d@(l eft) show the measured cross sections as a function of rriA 1 for the e + e 7 
and /u + /r “7 modes, respectively. The results include all data recorded at the Y(2S), T(35) and 
Y(4S) except for 5% data used to optimize the selection criteria. We estimate the significance 
for each fit by Ss = \/2 C/Cq where £ is the likelihood for signal plus background while £0 is 
that for pure background. Figures [5] (right) and [ 6 ] (right) show the resulting distributions for 
both modes based on 5704 and 5370 mass hypotheses, respectively. The largest deviation in 
the e + e - invariant-mass spectrum occurs at m ee = 7.2 GeV/c 2 with a significance of Sg = 3.4. 
Including a trial factor determined from the MC, the significance is reduced to 0.6cr. The largest 
deviation in the invariant-mass spectrum occurs at rn, J4l = 6.09 GeV/c 2 with a significance 

of Ss = 2.9. The application of a trial factor reduces this to O.lcr. 
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Figure 7. Exclusion region of the mixing parameter e as a function of the dark photon mass 
for present measurements. 


We combine the cross section measurements into a single result. Since we observe no signal, 
we determine Bayesian upper limits on the cross section at 90% CL assuming a flat prior. We 
include systematic uncertainties from the uncertainties of the dark photon branching fractions, 
luminosity and limited MC statistics. The resulting cross section upper limits are typically of 
the order of 0(1 - 10 fb). We then translate these results into 90% CL Bayesian upper limits 


on the mixing parameter as a function of the dark photon mass. Figure [7 


comparison to dark photon searches of other experiments: KLOE 49,50 
and </>(1020) —>• r]e + e~, electron-nucleon fixed target experiments APEX 
A' —> e + e _ , 

NA48 


shows our results in 
in e + e~ —>• /x + /r ~7 


the proton fixed-target experiment HADES 53 in A' 


51 


e 1 e 


and A 1 
WASA 


55 in 7 r u —» e + e 7 and electron beam dump experiments 56-58 . In addition, we show 
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Figure 8 . Expected exclusion region on e from results of new dedicated experiments APEX, 
DarkLight and HPS, MESA, VEPP3 and Belle II. 


the constraints from the measurement of the anomalous electron magnetic moment [59], The 
band of favored e versus taa' is obtained if the discrepancy between the measurement and the 
SM calculation of the muon magnetic moment 60 is attributed to dark photon production. In 


the 0.03 - 10 GeV/c 2 mass region, our results push the exclusion region substantially lower and 
supersede the results of our previous analysis 61 . Together with the results from NA48, we 


exclude the entire region favored by the dark-photon scenario for the g-2 measurement^ By 
searching for A 1 —> 7r + 7r , we can further probe the region near the p mass. Figure [8] shows the 
exclusion region of e as a function of niA 1 that is expected from future measurements such as Belle 
for 50 ab _1 and several new dedicated experiments including APEX, DarkLight and HPS 

MESA at MAMI in Mainz 65 and VEPP3 at Novosibirsk 66 


62 


II 

at Jefferson Laboratory |63,64 

Above 200 MeV/c 2 , the Belle II data will reduce the exclusion limit on e to 
lower masses HPS will push the limit below 10 -4 . 


2 x 10 4 while at 


3. Search for Long-Lived Particles 

The recent anomalous astrophysical observations have generated also interest in low-mass, long- 
lived hidden-sector states. Various models discuss this scenario, including those with dark 
photons (39,40,67,68 , an inflaton [69], supersymmetry [70j[7l], a dark Higgs boson [721 and 
with other dark-sector states [73]. Sever al multipu rpose experiments hav e condu cted searches 

and the multi-GeV/c 2 


56 


74 75 


76 


81 


mass regions. 


83 


for long-lived particles in the sub GeV /c 2 

In addition, dedicated experiments have been proposed [82] and some are under construction 
However, the 0( GeV/c 2 ) mass region has remained mostly unexplored. 

At the B factories, a hidden-sector scalar particle may be produced in an T radiative decay 
or in a B penguin decay and, in turn, may decay into a pair of fermions as shown in Fig. [9] So 


This assumes that decays to invisible particles are small, otherwise the bounds are weakened. 
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Figure 9. Lowest-order diagrams for the production of a long-lived dark scalar particle in 
radiative T decays (left) and in penguin decays of B mesons (right). 


we performed the first search for a long-lived particle L in the process e+e - —> LX where X is 
any set of particles and L decays to six different final states. Except for 20 fb _1 taken at the 
Y(4£), we use the entire BABAR data collected at the T(45), 40 MeV below the T(4£) peak, at 
the T(25), and T(35) corresponding to a luminosity of 489.1 fb -1 . We present the results in two 
ways. First without any assumptions on the production mechanism, we use the complete data 
sample and present model-independent results on the product of cross section, branching fraction 
and reconstruction efficiency of each two-body final state f, (jf = a(e + e~ —> LX) ■ B(L —» /)•£/• 
We have produced tables of the reconstruction efficiency as a function the L mass, transverse 
momentum pp in the CM frame and the proper decay distance cr assuming L to be a scalar 
particle. Second, we present results for B decays via B —» LX S where X s is a hadronic state 


with strangeness S=-l. This production mechanism is motivated by the Higgs portal 69||72 
and axion-portal [|84j models of dark matter. In this case, we report model-dependent results 
on the product branching fraction Blj = B(B —> LX S ) ■ B(L —> /). 

We use MC simulations of signal and background events to determine the signal efficiency 
and L mass resolution cr m . We generate signal events with EvtGen 185 and produce two types 
of MC samples. In the model-independent approach, we generate e+e - — > Lnir and B —> Lnir 
where n < 3. We generate 11 L mass values (0.5, 1, 2, 3, 4, 5, 6, 7, 8, 9, 9.5 GeV/c 2 ) assuming 
spin zero. This yields broad L-momentum spectra requiring efficiency tables to test specific 
models 


86 . For the model-dependent approach, we produce e+e 


BB, B —> X S L events 

where X s is composed of 10% K, 25% K* and 65% AT* (1680) [87]. The high-mass tail of the X s 
mass spectrum is suppressed by phase space limitations of heavy L-states. This choice yields an 
L-momentum spectrum as a function of m that reproduces the dimuon invariant-mass spectrum 
for B —> X s p + p~ |[851. In this approach, we generate seven mass values (0.5, 1, 2, 3, 3.5, 4, 
4.5 GeV/c 2 ). In addition, we generate the background processes e+e - —> BB (EvtGen 


85 


T + T~,p + p- (KKF2 [43]), e+e" (BHWIDE [42]) and qq (JETSET 
The detector response is simulated with GEANT4 [46]. 

We reconstruct L from two oppositely-charged particles (e + e~,p + p~ 
and J\+7r+) that originate from a common vertex separated from the IP by more than 3<r py 
where opy is primary vertex resolution. We reject background from Kg and A decays 
as well as low-mass peaking structures in the background by imposing mass thresholds: 


) where (q = u,d,s,c). 


, e+/T+ , 7T+7T , K + K 


me. 


> 0.44 GeV/c 2 , m efl > 0.48 GeV/c 2 


m h 


> 0.50 GeV/c 2 and < 0.37 GeV/c 


rn n7T > 0.86 GeV/c , mxn > 1-05 GeV/c 2 and mxK > 1-35 GeV/c 2 . The efficiency varies from 
47% for m = 1 GeV/c 2 , 2 < pp < 3 GeV/c 2 and cr = 3 cm to nearly zero for large m and large 
cr values. 

We extract the signal yield for each final state as a function of m using unbinned extended 
maximum likelihood fits of the m distributions. The strategy is based on the fact that a signal 









would produce a peak in the m distribution whereas background varies smoothly. Thus, the 
signal PDF is the pull histogram Pg{m) = Hg((mt — m)/a m ) produced from masses closest to 
the true mass rrif . measured mass m and its uncertainty a m . The background PDF P B {m ) is 
a second-order polynomial spline with knots located at the bin boundaries and is determined 
from the data. The bin width is variable and is defined as 15 times the RMS width of the signal 
distribution in that region. The factor 15 is chosen since it is sufficiently large to prevent P B {m) 
conforming to signal peaks but is sufficiently small to produce signal peaks from fluctuations. 
We scan the data in steps of 2 MeV/c 2 using the PDF ns ■ Ps(m) + ng ■ Pb^) where ns (tib) 
is the signal (background) yield determined from the fit. Figure [l0| shows the observed number 
of events in comparison to the expected background for each final state. Two mass points in 
the dimuon final state show yields with significance Ss > 3. At mt = 0.212 MeV/c 2 , we see 
13 events with Ss = 4.7. However, most of the vertices lie inside the detector material and the 
muon tracks have low momenta that are poorly separated from e and 7r. At mt = 1-24 GeV/c 2 , 
we see ten events with Ss = 4.2. We determine a probability of P = 0.04%(0.8%) of seeing yields 
with Ss > 4.7(55 > 4.2) in the entire dimuon mass region. Including the other final states, the 
probabilities are reduced by a factor of six. Thus, the data are consistent with background 
expectations. 



m [GeV/c 2 ] 


Figure 10. Observed mass distribution for T(45)+ off-resonance data (red solid points), for 
T(35) + T(2S) data (blue open squares) as a function of mass m for each final state with the 
background PDF Pb overlaid (red, blue solid lines). 
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For each scan fit, we determine the systematic error on the signal yield. The dominant 
systematic error comes from the background PDF. Other contributions are due to the mass 
resolution, luminosity, particle identification and the number of MC signal events. After adding 
all contributions in quadrature, we include the total systematic error (cr svs ) in the fit by 


convolving the likelihood function with a Gaussian function with width a sys . We determine 
uniform-prior Bayesian upper limits at 90% CL on df as a function of m for each final state. 
Figure [TTj (left) shows the results on 07 separately for the T(45) and T(35) + T( 2 S') data. Using 
the efficiency tables | 86 ], these limits can be applied to any production model. For the model- 
dependent interpretation, we set Bayesian upper limits at 90% CL on the product of branching 
fractions Bpf. Figure [TTj (right) shows the results on Bpf as a function of m for different values 


of cr. These limits exclude a significant region of the parameter space of the inflaton model 69 
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Figure 11. Upper limits on the cross section at as a function of L mass for each final state 
(left) for T(45) data (lower red curves) and for T(35 + Y( 2 S) data (upper blue curves) and on 
the product branching fraction B[ y f (right) for different decay lengths. 


4. Search for New 7T° -like Particles 

The measurement of the 7 r° form factor in two-photon collisions by the BABAR Collaboration 

At sufficiently high squared momentum transfer Q 2 , 


90 


93 


89 


the 


94 


has raised various discussions 

pion form factor should approach the Brodsky-Lepage limit of V^fn/Q 2 — 185 MeV/Q 
At Q 2 > 15 GeV 2 /c 2 , the data is expected to be described well by perturbative QCD. Our data, 
however, show no sign of convergence towards the Brodsky-Lepage limit as Fig. 12 indicates. 
Though the Belle data 95 show a better agreement with the perturbative prediction, they are 
consistent with our results. So, maybe a new particle </> exists with mass close to that of the 7 r° 
decaying to 77 and thus causing the non-convergence of the ir° form factor. The new particle <f> 

or a so-called hard-core pion (n® hc ) 


96 


may be a scalar ( (j>s ), pseudoscalar (4>p) 

We focus on a search for 7 r°-like particles that are produced in association with a t + t~ 
pair 1971. Figure [l3] shows lowest-order Feynman diagrams for the production of 7 r°-like particles 
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Figure 12. Measurements of the 7r° form factor as a function of Q 2 . 


in association with t + t . For Q 2 > 8 GeV 2 /c 2 , predicted cross sections for e + e —> t + t 4> are 
large: cr n o = 0.25 pb, a^p = 2.5 pb, and = 68 pb. Thus, we expect large event yields in 
our data set: ^ t 7T ° hc = 1-2 X 10 5 ,1V^ P = 1.2 x 10 6 and N ( j >s = 3.2 x 10', respectively. 


i 

n° k 



e r 



Figure 13. Leading-order SM diagram for t + t annihilation into a 7r°-like particle (left) and 
for the radiation of a cj) from a r lepton in e + e~ — > t + t~ (right). 


Using the full BABAR T(45') data set of 468 fb -1 (4.3 x 10 8 r pairs), we select t + t~ events in 
the decays r + —> e + u e i/ T and 1^] decays requiring px > 0.3 GeV/c for each lepton. 

In addition, we require exactly one 7r° —>• 77 with energy 2.2 < E n 0 < 4.7 GeV in the laboratory 
frame. After excluding the energy of the photons from the 7 r° decay, we require the remaining 
energy E extra in the electromagnetic calorimeter to be less than 0.3 GeV. We reduce background 
from radiative Bhabhas by imposing a minimum photon energy of > 0.25 GeV and an 
opening angle between e^ and 7 of 30° < #(e, 7) < 150°. We further reduce background from 
semi-leptonic and hadronic decays by requiring E sma u + E^ > 0.5 Eqm and m n ± n o > m T 

where E sma u is the energy of the lower-energy track and Eqm is the CM energy. Figure [14] shows 
the Es ma ii + Eff, distribution before the requirement of E sma u > 0.5 Ecm, which is consistent 
with the expected e + e _ —> t + t~ energy spectrum. We simulate background modes e + e _ —> BB 
with EvtGen 85 , e + e _ —> qq continuum with JETSET 
with KK2F 


e 1 e 


R R ,T'T 


88 


e 1 e 


e 1 e 


with BHWIDE 42 


43 where we use the TAUOLA library 98 to generate r 
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charge conjugation is implied unless stated otherwise. 
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decays. We model radiative corrections with PHOTOS 
response with GEANT4 [461. 


99 and the simulation of the detector 



Figure 14. Top: the E sma u + E$ distribution 
for events with 100 < m 77 < 160 MeV/c 2 in 
data (points with error bars) and simulations 
for e + e _ — > t + t~ , /r + /x _ , bb, cc (colored 
histograms). The dashed histogram shows the 
expectation for i t q hc production at a = 0.254 pb. 
Bottom: pull distribution. 


points with error bars), expected sig¬ 
nal (dashed histogram) and simulated 
e + W —> t + t~, bb, cc background 

spectra (colored histograms). The inset 
shows the low-mass region in bins of 
10 MeV/c 2 . 


Figure [15] shows the entire m 77 invariant-mass spectrum after applying all selection criteria. 


For comparison, we also show a simulated signal and the simulated spectra for e + e _ —> 
t + t~, bb and cc backgrounds. We fit the m 77 invariant-mass distribution in the 

region 50 < m 77 < 300 MeV/c 2 with a Gaussian signal and a linear background: N(m 77 ) = 


N, 


m. 


' 77 1 


We use extended unbinned maximum log-likelihood fits 

2 


g G(flg, (7 g) T IV/(l + 

to extract Nu n ,N g and au n . 
determined from control samples. The mean jig represents the mass of the <j) particle and is fixed 
in each fit. We scan n g values between 110 MeV/c 2 and 160 MeV/c 2 in steps of 0.5 MeV/c 2 


The width a g of the Gaussian function is fixed to 11.1 MeV/c 


and extract N„ from each fit. Figure 16 shows the resulting N„ distribution. Figure 17 shows 


the fit for the highest yield. We observe a raw yield of N g = 6.2 ± 2.7 ± 0.06 events at a mass 
of 137 MeV/c 2 . After subtracting 1.24 ± 0.37 events of peaking background and correcting 
for a fit bias of —0.06 ± 0.02 events, we find 5.0 ± 2.7 ± 0.4 signal events. We determine 
identical efficiencies of (0.455 ± 0.019)% for tt q hc and (j>p. For ^> 5 , the efficiency is reduced to 
(0.09±0.004)%. Systematic uncertainties include contributions from the MC sample size, energy 
scale, 7T° efficiency, energy resolution, particle identification, momentum scale and momentum 
resolution and amount to 4.2% for tt q hc , <j)p and 4.4% for This leads to measured cross 
sections of cr „o ^ = (38 ± 21 sta t ± 3 sys ) fb and = (190 ± 100 sta t ± 20 sys ) fb, respectively. 
Since the result is consistent with background, we set Bayesian upper limits on the cross sections 
assuming a flat prior. We obtain <r^ H t< j> p < 73 fb and < 370 fb at 90% CL. This should be 
compared to predictions of 270 - 820 fb and 68 - 1850 fb, respectively. 

To check the compatibility of the measured production cross sections with the contribution of 
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particle (f> to the 7r° form factor, we include this measurement as an additional term in the x' 2 of 
the fit to the Q 2 dependence of the 7r° in our data. The increase in \' 2 follows a x 2 distribution 
with one degree of freedom. Thus, we can determine p— values for the process e + e _ —> T + T~(j) 
to contribute to the 7 r° form factor. We obtain P 7T o hc = 5.9 x 10 -4 for Pcf> p = 8.8 x 10~ 10 

for (j)p and p^g = 2.2 x 10^ 9 for efts- We conclude that a new particle 0 with mass close to the 
7T° is not a likely explanation of the excess of the 7r° form factor seen in our data at larger Q 2 . 
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Figure 16. Number of events N g as a 
function of the </> mass hypothesis p g 



Figure 17. Measured m 77 spectrum (points 
with error bars) with fit overlaid (solid line), 
peaking background contribution (dotted line) 
and expected backgrounds from e + e _ —> 
t + t ~, /r + /x _ , bb, cc (colored histograms). 


5. Conclusion and Outlook 

The B factories are an excellent laboratory to search for light dark matter. We see no dark 
photons in the 0.02 — 10.2 GeV/c 2 mass region pushing the limit on the mixing parameter e 
at a level of 10~ 4 to 10 -3 depending on the dark photon mass. We performed the first search 
for long-lived particles at a high-luminosity e + e _ collider and the first search in a heavy-flavor 
environment. We observe no long-lived particle in the 0.2 < m < 10 GeV/c 2 mass range and 
proper decay lengths of 0.5 < cr < 100 cm. In a model-independent approach, we set 90% 
CL upper limits on the product a(e + e~ —»• LX) ■ B(L —>• /) • e/ for different masses, pp values 
and proper decay lengths. We also give model-dependent 90% CL upper limits on the product 
branching fraction B(B -x LX S ) • B(L —> /). Finally, we searched for the production of 7r°-like 
particles in association with a t + t~ pair and set upper limits on cr(e + e~ —> t + t~ 0) < 73 fb 
at 90% CL. The hypothesis that a 7r°-like particle causes the non convergence of the pion form 
factor has a p-value of 5.9 x 10 -4 for tt q hc and thus is rather unlikely. For cj)p and 4>s production, 
the p- value is even lower. BABAR will continue with searches for new physics. Belle just released 

• A'h', h! 


first limits on e + e 


A'A' 


100 


Since Belle’s luminosity is about a factor of 
two larger, the limits are a factor of two lower than the BABAR results. However, substantial 
improvements are expected from Belle II. At 50 ab~ 4 , present BABAR limits will be improved 
by up to two orders of magnitude depending on the final state. 
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